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ABSTRACT 


A threshold detection system based on the Neyman-Pearson 
criterion is derived for an infrared nutating optical system. 
Detection is optimum for small signal-to-noise ratios and) 4 
Gaussian uncertainty in the pointing error Ol ties ope lea! 
System. The signal spectrum and background power spectral 
density for a rectangular space filter are computed numeri- 
cally and used to specify the matched filter for the thresh- 


old detection system and evaluate its performance. 
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I, INTRODUCTION 


Pi itsenared (IlRj@deteetien systemedetermines the 
Meesenice Or absence of a target in a particular area of 
mee by processing the received IR energy from that area. 
The target must be found amid ever-present background 
radiance and ay Internal Nomse Peneratedubywethe detection 
meee. In enis thesis, one type of system, called a 
Meeating optical system, is considered. 

iMirared detectiom 1S a Combination Of Spatial frequency 
processing and temporal processing. Spatial frequency process- 
ing 1S a means of achieving a modulated output from the IR 
memance passing into the optical system so that information 
ieee Casily extracted [Ret. 14]. To extract this informa- 
Mmiom., some form of temporal processing) is used. In a 
Memeatane optical system, the optical axis rotates about an 
axis perpendicular to the image plane in such a way that a 
POoMlemiinage traces out a circle in the image plane. The 
Spatial frequency filter ioe rCeerO tain esensdel Ve Mace m ia E 
tiene image plane which produces a voltage output that is 
@meiitied and fed into a temporal filter for processing. 
Pemutating optical system is frequently Used in a missile 
feexcr because it produces a signal when there is no 
@raeking error and the rotating optics can be made into a 
gyro that provides the inertial reference. 

This work is a synthesis and extension of two recently 


published articles in the literature. The first, by 





Heegcer ket. 3), formulates the detection theory of a 
leew signal in background and white noise. Since in IR 
@emeerecdetection the signal is never known completely, 
bom S wWOrk 1s extended to the case where the target 
amplitude and position are unknown. However, a known 
Mueeet Shape 1s assumed as well as a priori statistics of 
Mimacude and position. Since detection 1S most difficult 
when the amplitude is small, a test (called the threshold 
detector) is derived that is optimum for weak signals. 

lites pCeiiaune storm ot the thnesnold detcetor and to 
meremlate 1(S performance, an eigenvalue integral equation 
must be solved. Ordinarily, this is a difficult undertaking 
Peeecme SOlution becomes trivial if the covariance function 
fmeme background process is periodic. This is exactly 
emenwcondition realized for the background process out of 
a nutating detector provided that the phase, which is 
physically unimportant, is averaged out. 

Piuticlssonm' s results (Ref. Iii specity the form or the 
Signal and background coefficients in terms of the system 
pomametecrs and signal and background radiance. These 
derivations were first published in a Swedish internal 
mepent (Ret, 10], then re-derived in English by Yoshitani 
[Ref. 6]. The latter's derivations are included in the 
appendices . 

samuelsson’s equations have been applied to a rectangu- 
lar IR detector assuming a Wiener spectrum for the random 


background and the coefficients have been calculated. 








lemmestiese CoOecti1cients, the optimum filter has been 
determined and the probability of detection computed as 
a function of signal-to-noise ratios and background-to- 


morse ratios. 





iene hei) ROCESS ING 


ites patialeunoqucneys:1lterineg portion of the detection 
system is comprised of a lens system and an IR sensitive 
meeector located in the image plane, the focal plane 
on the Opies, io adequately derive information 
fmeomecne IR radiation, the lens system and detector must 
be of a nature that causes the voltage output to be modulated 
memeeway that information about a target within the field 
@eevtew can be processed. In a practical nutating optical 
eerem, the lens system moves and the detector 1s stationary. 
Pewever, for conceptual purposes, this 1s equivalent to a 


Beatdenary Lens systen and nutating detector, 


A. TARGET SPECTRUM 
Tiemoomeet tOuberdetecetca semertarcet. 1S located 


Marne Object plane a distance R from the lens and parallel 
to the image plane as shown in Figure 1.1. The angle from 
micmepenpendicular axis connecting the planes to the target 
Soomimates. measured in the x-direction is called x and is 
measured in milliradians; the angle in the y-direction is 
femeuire target 15 considered to be close to the optical 
peota and at a Sufficient distance, R, such that small angle 


approximations may be made, 
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Consider a target with a radiance distribution 
N(x) (W/cm*-sr) HEC mrss tne —cocanalt power into a Unitysolid 
angle per unit projected area of the sources. r mes) 1S 
MevectOr whose origin is the perpendicular axis connecting 
emer coordinates. 

Assuming no atmospheric attenuation (this can be 
mmemuded in N(r) if it exists) and perfect optics, the power 


Oiler ibution in the image plane is 


P. (r) = N(r) * Ag (w/sr) C229) 


where Ag tSUeheweCtbkeetIVe Chtmance areameot the opties. 


Divadimes by A one can see that the radiance 


Orc 
iietribution in the image plane 1s numerically equal to the 
maemamee diStripution in the ooject plane. 


HeWevelmweDecause tiemoprnes 15 NOt pPperiech, a 


peat Object given by 


Hp 6(r) = Hp 6(x) é(y) Ww) , (2. ) 


(where Hp (W/em*) 1S the Lrrddiance arty the GOptrtcs Trecel ved 


Mmeone tine point source) iS imaged as a power distribution 
Hp + £0 (1) (W/em’ Sr) cn> 


where f,(r) (sr 735 imine pOm ims tread shunceTon Of Lie 
Bertcs. Irradiance is defined as the radiant flux indicent 
fommer SUTrface Of unit area [Ref. 5]. 

themenmore, thewradiance distribution, NiGor in, tne 


object plane will be imaged as 


N' (r) = [N(r-s) f)(s) d?s (W/cm*Sr) (2.5) 


ee 
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where dee doa sic uineminneOrat1On 1S over the entire 
image plane. 

thesdetcetor scea Space filter in the image plane 
which may be a moving retical or simply an aperture across 
waren the power 1S integrated. 

The power in the image plane incident on the detec- 


ion 1S 
H = [N'(r) t(r) d’z (W/cm*) (2.6) 


Wteac 1(r) = t{X,y) 1s the transmittance of the detector. 
When the detector moves as a function of time, the 


Pemer incident on the détector becomes a function of time, 
z | Z eee 
H(t) = [N'(r) t(r-e(t)) dér (W/cm?) ae 


fimeeer ott) describes the movement of the detector in the 
iiewee plane. 

PeViut Ine GebCctOr MOVesmwcIneular ly an the image 
meme DUt 1ts Orientation remains fixed, i.e., the coordinate 
eruem) (X' ,y') Shown in Fig. 1.2 does not rotate. The nuta- 
Meoneradius measured with respect to the optical axis is 
given by o. 

Assuming a nutating optics and no relative motion 
Peewecen the target and the optics, the radiance on the 
WemectOr will be pertodic and thus can be expanded in a 


RPourier series 


H(t) = ZH e Chee ee (2.8) 
ial 





where 


T -jnwot 
i H(t) e Git 
0 


ao 
M 
sel oo 


Pemlacing H(t) by equation (2.7), 


1 a -jJnwt , 
Be = Lak dt e d- maior)... (r=o(t)) (622-9) 
i ie is iN 
ihe motion deseribine function, ole LOr weir ou lar 


fer ation 15S 


O Ci )its= p(cosw.t, Sin wot) ; (25 18) 


where Has 2u/T (rad/sec) is the radian frequency of nutation 
aoe t 1s the period for one nutation. 
SUbSt2butine (2.10 )81n (2.96 Me can be shown (see 


Appendix A) to be 


th, *) a? 


H = 5" fuck) r*(k) e279 J (2m0 \ k (2.11) 


CiyCibar?) NUtat 1OnM oa 


— LO y' Sat ial ee 







ak —— GueiiGy sat ter 
a >» (detector) 
V7 \ 
/ : eens 
IMAGE OTe aneee ye 
; Space filter_” 
/ ae coordinates \ 
| oe oe 6. = Nutat ion mead. 
Image plane p 
coordinates 
X 
Preure dye.) CireulareNuGation in ithe Image 2 lanc. 


14 





where k = (k 


k_)}, two-dimensional spacial Fourier 
x? y ’ £ 


transform coordinates [Ref. 8], 
N'(k) = transform of radiance diiSserupueion in the 


image plane, 


mK ) complex conjugate of the tromesormmeot the 
optical transmittance function iit themdeteetror 
@oOa quia cesSyStcin, 


a aul? 
@ = tan KP Ky: 


Mueate ron radius, 


o 
i 


and J (2) jottheemtm order Bessel function of the rirst 
kind. It has been assumed that the image plane is very 
large compared to the radiance distribution from a target 
so that the image plane may be considered infinite for 


mathematical purposes. 


B. CORRELATION OF BACKGROUND NOISE 

The background power incident on the detector is from 
a sample background scene where a scene is a two-dimensional 
random process Breteteriaen by the radius vector r = Coy) 
from the optical axis. 

pee N(x) be a sample scene radiance distribution on the 
object plane from an ensemble of scenes which has been given a 


emlcable probability structure. The background power incident 


on the detector is 
B(t) = J N'(r) t(r-p(t)) d’r . Oa 


The covariance function is defined by taking the expected 


value 


is 





Pete et) } 18 (oO) El Ba) )] = Rp (tu) (2513) 


where E denotes the expectation over the ensemble of scenes. 
Meeause the detector is nutating, the random process, B(t), 
will not be stationary. However, if the detector and 
background do not move with respect to each other, then 


Rp (t,,t) Po VAdOUpD I Petr Lod le. 


Pc tite Ni) = R(t 5c.) (2.14) 


where M and N are integers. 
It can be shown [Ref. 7] that the form of the covari- 
mice Of a doubly period process is 
00 Jmw,t 


Rp (t},t>) = FY 86 e 


{| a INWet?2 
mn 
m,n 


(25 ES) 
Because the nutation 1S periodic, over the ensemble of 
Scenes, there is little significance in the starting time t. 
Therefore, considering t as a random variable that is 
imetorm over the nutation period, one Gan “average it out” 


piece cOVariance function as 


T 
p(t) : - Ra(t+t,t) dt . aakee 


Pliestituting (2.16) in (2.15), 


00 jmu , T J(m-njoot 
ae ait (eelie 
0 


! 
ee | 
WD 

= 

2, 

o 
| 


Rp (7) 
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feane the relation 


T J (u-vjayt 
ac cl eer [ 2:ehSs) 
0 


= 


6uv = 


f@iae correlation function becomes 


© JNW yt 
Ry (v1) = : BL eC ; C2) 


Ea@uation (2.19) shows that B({t}) can be considered a 
wecie) Sense stationary process and that it is periodic in 


meyemmean Square sense; 1.¢., 


Ry (1) = R, (T+T) for any tT (2220) 


iMmiespower Spectral density 1s the temporal Fourier 


transform of Ry (vt). Taking the transform one has 
ee = 27 * BO 6 (w-nwy). (2221) 


The coefficients, 3,» can be related to the optical 


System parameters by 


2 ee” 


2 es 2 2 Z 
8, = flsGl® [Fy (KL Wa (kK) J° (2me"K +k") do (2.22) 


where Fy (k) Sune  transctorm 01 the point spread. tumeti on, 
and Wp (kK) is the Wiener spectrum, the transform of the back- 
@mpound Correlation function. This derivation was first made 
by Samuelsson [Ref. 10] in an internal Swedish report then 
again by Yoshitani [Ref. 6]. Yoshanti's derivation has 


been included as Appendix B. 
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bogatton (2,19) implies that the background, B(t), can 
be represented by a Fourier series with uncorrelated 
@eerticrents [Ref. 7]. Thus, 
JNuyt 


B(t) = : bY e (2a2s) 


where the convergence is in the mean square scnse. 


ie coefficients, ba? Satisfy 


ie (Baa) n= 0 
EDe) = (224 } 
n # 0 
Memeover, they are uncorrelated 
a m=n 
Co 8 = : (225) 
0 m7#n 


een llTE NOISE COMPONENT 

Mian actual system there will be a certain amount of 
Mnaeommal molse generated in the detector and preamplifier 
which is additive to any background noise or signal. In 
feomrreducencles OL interest the primary source of this 
memoe 1S thermal agitation. This type of noise, called 


Momrson noise, is assumed to be white and Gaussian. 


Let W(t} be Gaussian white noise with power spectral 
density ‘I (V*/H,). The power spectrum at the output of the 


preamp with bandwidth B is shown in Figure 1.3. 
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on 


PE hoor OWwen Spee CMiuniwon int te Nose . 


iiie RMS voltage measured at the™output of the preamp is 


bof 2 


Vrms - - 2B = NB (v*) : (2220) 


Is, 





i ee eRe hae ROCCE Se N.G 


Iiemoutput voltage from the detéetor is temporally pro- 
Serecdweo determine the presence of a target. The type of 
pmecessOr 1s determined by applying statistical detection 
theory to the known behavior of the nutating system and var- 
memeerealistic assumptions about the system and the target. 
mop timum processor in the sense of small signal amplitudes, 
eed a threshold detector, is developed in this section. 

Pomerer weet. 5) derived the detection theory for the 
how Signal case in the presence of additive background and 
Moeeeeenolse components. In IR detection, the target ampli- 
midewand position are usually not known completely. The 
threshold detector is an extension of Harger's work which 
feecountsS for the unknown amplitude and position. 

ie form of the optimal detector under the Neyman-Pearson 
Criterion will be derived for multiple observations. It is 
assumed that the image plane and object plane remain parallel and 
Seoeronary with respect to each other during the observations. 

The detection of the target under the Neyman-Pearson 
criterion becomes a problem of hypothesis testing. The 
task of the detector is one of choosing between two 


hypotheses, H, that only noise is present and Hy that in 


0 
addition to the noise there is a target present. The design 
@eethe detector is one that permits the correct choice of 
hypothesis Hy (a detection) with maximum probability of 


detection, Qd, in a fixed probability, Qa: of choosing 


20 





Hy when Ho 1S true {a false alarm). The structure of the 
Meeman-Pearson criterion requires forming a likelihood ratio, 
Ajo» which is compared to a known threshold, Ve Specified 
fiemeea given false alarm probability. If the likelihood 

Tatio exceeds the threshold, hypothesis Hy 1S assumed true 
and if not, hypothesis Ho iswassumed@Gmic. The logic process 


may be written 


H 
> 
10 < V can 


mnere the threshold, Ve is derived from the expression 


©O 


Rea = f p(Z|H,) dz (Sez) 
L 


mmmomis a SUtticient Statistic of the received data. 


A. DETECTION OF A KNOWN SIGNAL 

Meare data fOr processing 1S a set Of Me functions of time 
oat), n=1,...,M} where each function is the output of the 
meeamp during one nutation (0,T). T is the nutation period. 


The output may be represented as 


Z(t) = S(t) ae ete) oe CE DE Se ols roel Cee) 


S,, (t) represents the output due to a target. The component, 
Pht), represents an additive background and is a sample 
Mmmetion ©f a random field (B). The relative size of the 
target in the field of view is assumed sufficiently small 


that the background is additive while B(t) is assumed to be 


Dak 





| | a, 





the sample function during each of the M nutations observed. 
at) MepLesentsS the additive nelse component due to the 
mmcernal noise of the detector-preamp combination. It is a 
sample function of a white Gaussian random process Cis of 
Zeno mean and power spectral density N/2(V°/H,). 

Mie two hypotheses for the Pikelihood ratio may be 


Written as 


H,: {Z(t ) S,(t ) + B(t ) + Ny (t yy Cade ee ahi 


(3.4) 


Hy: Ge) seas: Ne me Cercle See 


Denote the likelihood ratio as Ay) ({Z,}) where the like- 
lihood ratio is defined as 


p({Z }/H,) 


NygMZ,3) = p({Z_}TIH 


(o5) 
n 0) : 


Mace likelihood ratio 1s complicated by the fact Zn and 
ie are not independent because of a common noise component, 
moe lO simplify the derivation Harger introduced an 


avo l lary hypothesis 


Hee ec = NO ymin ann, (3.6) 


oe rel 


mnen the likelihood ratio, Ayg> Can be computed using 


jee chain rule for likelihood ratios: 


eee v eC TM CAE Iii) Ny oe) Carts nh) 


"10 ~ pClZ Fy)  PCZ_FH))/PCZ,} | Ay) 


= Ayo/Koo- 


Pee 


If the sample functions for the background, B, were non- 
randon, Z.(t) and Z(t), i # j}, would be mutually indepen- 
dent because of white noise. Accordingly, if B were fixed, 
@me can write the conditional likelihood for the ae 
observation Mie Crean rr Moreover, the likelihood ratio 
fOr M observations would be the product of the likelihood 


menos for €ach single observation, 
: n 
ea Paz SER IB 


The conditioning is then removed by averaging over the 


ensemble of backgrounds 
Ayo (CZy}) = By[My 7 (tZ,11B) | 


Eemree the noise is white Gaussian, the conditional 
Heeellihood ratio is that for detecting a known signal 
Se, in white Gaussian noise and has the well known form 


eet. 15), page 253) 


M cy 
yy ee Z(t) (S_ (t) se sbulie lhe toe 
n=l 0 


ZN 


Aya: ({Z,3|B) = exp | 


Z |e 


T 2 
Gre Aron ac] 
0 


men may be rewritten as 


M af; 
Nn a 2 oy ¢ 
Ay> (Z| B) exo | 3 — ; Z(t) S(t) at 
1 M T 9 
-— y ff Ss, (t) dt} exp a(b) (Cae 
N n=1 0 








where 


a(b) = 


ius 


Pe 7 Gen es Pie ycues ee par. a 
f (Z(t) + S,(0)) Beat - FL BA Cr)ae, (3.15) 


n=l] 


. n 
To find A, )(Z,)), one must evaluate E,[A 12 (12,3 1B) J 


where 


T 


M 
E,{A",,({Z,}|B)] = exp {f zs 2,(t) S,(t) at 


pl 
qh 


ii 
(3. <1e6,) 


T 
: s “(t)dt] - E, Ges ob). 


Z| 
Te cee 
b-4 


Assuming the background is Gaussian the expectation can be 
evaluated by representing the background in a Karhunen- 


Loeve (k-L) expansion 


B(t) = : by aes. ee Ua basal) 


where >, (t) boweiMe cle enTunctiOn Commeesmemarne (fO Cicemy asic 


dy of the homogeneous integral equation 
| iG 
Ayo, (t) = ; Rp (t,u) >, Cu) du wel eee Sas) 


Rp (tu) is the covariance function of background 


Rp (t,u) =bd (et ie :- Mp (t) } [B(u) - Mp (u) J} 3} (3.19) 


where Mp (t) is the mean of the background. 
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aiinte 1a En lexpa(b)) can be calculated in a straightforward 


manner. After the lengthy calculation shown in Appendix C 
In E, texp a(b)} = - i ewe CL az ze rA,) 
B P 2 Te 
k 0 
> iT M T M 
+ N eat z (Z(t) = Ss (t)) f du : (Z_(uj)-s,,(u))q(t,u) 
0 n=l 0 m=] 
1 it 5 M 
- 5 : ements) ne, (Z_(t)-S, (t))} h(t) (3.20) 


where q(t,u) satisfies the integral equation 


ih 
BS meu) + {J Ri(t,t)a(t,ujdt = R,(t,u) Or ue, <P Zn) 
2M 0 B B — — 


acim mit) satisfies 


z 
ime) + Of Ry (t,t)h(t)dt = mp (t) 0 eS ee er (S22) 
0 


=— 


2M 


Sees) is the impulse response of the filter that gives 

Heeeminimunm mean square error estimate B of the background 
Mien the mean is zero [Ref. 15]. The additive white noise 
is reduced in amplitude by 1/M. Using the K-L representa- 


meron, one can also write 


. 


Ax 


Le) 
—_ 
— 





a(t,u) ~ 6,(t) ¢#(@u) O<u,t<T = (3.23) 


ae, 


aa Sad uae 
Sz 


k 


Pewee wequat ions (2.9),  (3:10)5-(5.14) “and (3.20) and 
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observing that 


i cee Mette Ds. =" 0 (3.24) 


it can be shown (see Appendix C) that the likelihood function 


pepenypotheses H, and H, is 


1 0 
M T 
= t 
A, 9 (12,4) exp { 2 if Z(t)q_(t)dt + Kt} (S250) 
n=l 0 
where 
9 1 T OM 
q(t) =35 [S,(t) - oy ip y S, (ujq(t,u)du] (55.25) 
O r=] 
and 
reesei M T (C , ITM M 
Meee 6) SS (tidt + ff EF We See Ge CER eGeiel))relretel 
n=10 % Ne ba 
, M T 
-— = & f SomGe ntact oe 3) 
no. n 
n=l 0 


ime 1OSdrithm 1S a monotonic function. thus it offer 
Momloss Of information but greatly simplifies the calcula- 


tions. In the remaining work 


A 


‘10 (123) = In Nig ({2 3) (35-20) 


Rewriting equation (5.25), 


At, ) (2,3) =e (123) ase (S25) 


where 
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it 
; Z (tia, (jdt 


ence ee 


Ht 


K' does not depend on the input data so it may be included 
in the threshold. The hypotheses test now becomes 
ay 
L( 2.) 2 Threshold . (3.28) 
Ho 
It is also possible, as Harger has shown, to derive 
tests similar to (3.28) without making the Gaussian back- 
Meena assumption. In particular, if the background fluctu- 


ma-emiiany times over (0, T) such that a(b) of (3.15) is 


Momoximately Gaussian, then the test is 


H 
Meet : 
> 
heel Z (tig, (tjdt . Threshold (5.29) 
n=l 0 H 
0 
where 
. it M 
sme 2M 1 
g,(t) = [S,(t) - SPS Ry(t.u) ff 2 S,(u)]dul 
0 Q=1 
Since by Mercer's theorem 
Rp (t, uv) = : Ayo, (t) o¢ (4), US te, Wea cael (S250) 


on comparing ge: &. ome (52) )elt 1S) ScoCn that ene rabove 


model approaches the Gaussian model if 


a 








iitenebove Condition, which implies a weak background, is 
sufficient but by no means necessary. The two models 
approximate each other whenever the second term of qn and 

g, are Nearly tacmiical, Ine point 1s that while the 
Gaussian background assumption may not fit the reality, 
mmpocesSor based on the assumption is probably not far 

from being the optimum. Also, the Gaussian assumption allows 


one to compute the processor performance. 


Bee DETECTION OF SIGNAL WITH UNKNOWN PARAMETERS 

Mies likelihood ratio derived in section A 1s applicable 
Cie eit the signal is known completely. Such is seldom the 
@ioenin IR target detection and in general the signal 1s a 
function of some unknown parameters. Typicaily, the ampli- 
fmig@ewa,Of the target and its position ie iene IMace lan 
are the unknown parameters, and the form of the signal can 


be rewritten to include these parameters: 


- 

> =: S(tia, Te) 
mir thermore, by assuming that only point targets are of 

interest and including the assumption of fixed image and 

object planes over the number of nutations of interest, 

the shape of the waveform will be the same during each 


Mutation. However, the amplitude is allowed to vary between 


emecesSive nutations. 
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The signal waveform may now be written as 


> > 
S(t 34,7) =o f(t,r,) (S254) 


where a iS) the sitemal amplitude and ea) 1S) the steal 
shape for a target at Tie 


Rewriting hypothesis H, under these assumptions 


i 


gives 


Se Te UGS) e@ueg) + B(t) + N (t); n=1,...,M 


OR ata ls (sae) 


ime likelihood ratio Mio becomes the conditional 
likelihood ratio Nyse ene) obtained by replacing 
S(t) by ese) everywhere in (3.25). 

Assuming a priori statistical knowledge of a and —, 
and their independence, the likelihood ratio may be 
averaged over the respective density functions to produce 


apeencONditional likelihood ratio 
A, ({Z,}) = [{moCizytla. top lalp(t,)dadr,. (G55) 


if e > 
Imstead of averaging as above, one could estimate a and ue 


by the maximum likelihood procedure and substitute the 


(aS 


ne > =. oe a >. = —— _-_ YF ~~ 





- rr en ae - > 7 


: =a e e e « e 
estimated values a and ry in the conditional likelihood 
m@aero and try to maximize the generalized likelihood ratio 

Re 
Ay QZ, t]a,r9). However, this procedure is not followed 
here because the maximum likelihood estimates are difficult 
hemobtain. 
Since the main interest 1s in detection of the weak 
mommeeet, the small signal case is of primary interest. 


Expanding Reyer2a fann} in powers of signal amplitudes 





M 
-> _ A) A -> 
Ay 2,3 a,TQ) = A ee on Ja Ag QZ tlasrg) lang 
. Seon) 
l M M 52 | | 
Te a? ee i et 2 ae 
Zee ie a Jada, FO no '~? Q*'a=0 


imcecratine (2.34) as in (3.33) term by term and 
Pertaining only the first non-zero term that depends on the 
data provides a statistic which is optimum for weak signals 
and is called a threshold detector [Ref. 4]. 


The evaluation is made easier by the fact 


_o 
Jay 





3 ae 
pemee 10627) 3% G)) = 


irom CE | aise) (S335) 
k a=0 ae 0 


az; 


which results from 


> 
A({Z 3 }0,%) =i 
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The first non-zero term becomes 


rm 


0 a 2 ; > 
mee A. ,(1Z.}l|a,r,) | Sal VL: tlieeeen: .oidlet 
day 10 gee wl eres, Nop on 0 
a T a 
Pa if Z (t)dt i q(t,ujf(u,rp)du 
0 0 
a Gace 
- f £(t,r,)h(t)dt 
0 


e Sd alg 
Pemeomnineg the imtepgration over r 


0 
a : i : 
A, ) (12,3) = i & v 7 ale Z(t) Ge) os qt eu) en Gurca 
n=l 0 0 
T. 
-f £(t) h(t) dt (3.37) 
0 


mieweeeune bar denotes the averaged quantities. The first 
Guemiast terms are independent of the data and may be 
included in the threshold value. Moreover, if the ampli- 
midiemarc identically distributed, as 1S generally the case, 
wicmmmreshotd detector performs the test 


H 


M oT a 
We = eee ce) Ct) ad & : W (3.38) 
noo u 
H 
0 
where We 1s the threshold and 
- Dee qT : 
q(t} = vy Lf (Ct) ea te Une ain) te 65 352) 
0 


Eme 





Thus the threshold detector can be mechanized either as a 
@Sourelator or matched filter followed by an accumulator. 
The filter function q(t) may be represented in a series 
eapamsion in terms of the eigenfunctions of Ry (t,u) and the 
Seerricients of the signal expansion. Expanding f(t) in 


terms of the eigenfunctions: 


Gas 2 eee) (3. 40) 


k 


eiesotttuting (3.40) in (3.39) and using (3.23) 


— 


ie 
- N 
i es ON 


rem (te) 5 


k 
+ Ay 


a 
2M 
C. EVALUATION OF PERFORMANCE 
The Gaussian assumption was made in deriving the 
Mepesnold statistic; therefore, the threshold statistic 
mescift 1S Gaussian. This means that tne probability of 


detection, Qd, and probability of false alarn, are 


Qra? 

laure ly determined by second-order statistics, e.g., 

the means and variances of W under hypotheses Ho and H,. 
Recall that the threshold statistic 1s 


je 


Te se 


n=1 


Mire smeans @nid Varrances ot W are: 


SZ 


o, () , (5) 


f Z(t) ait at: (Se) 
0 





M al 
E(W[H)) eee wece Eee) + No (t)) dt 
me Nl 
Se 
M TT. 
= »y f q(t) m.(t) dt 
n=l] 0 B 
MoT. ‘ 
B(WH, ) = eeeecie ye emer (t.r-) + Bit) + Ni (t) ] dt 
és n 0 m 
n=l1 0 
M  T_ ; 
== E(WIH)) ele wert leatt) dt 
n=1 0 


VAR (W[Hy) = VAR (WIHT, ) = VAR 


nm - -2 oe - 
=~, f q°(t) dt + M’ sf q(t) Sf q(u) R,(t,u) dudt 
| ee 0 0 2 

(3.43) 


Letting p,, (1H, ) be the Gaussian density function of 


W under H. Vineieie 1S SpeCd (le aban co mC ell au lola db OMe. Nuilcnl 


= PP Coie (3.44) 
W. 


QrA 


and 


Oh) =f penesiine ).dx (3.45) 
les 


where ie is the threshold voltage. 


E@iaetonS (9405) mand (4,44) can be solved to yield 


a 


a 





Qa =O. (x) (3.46) 
and 


Qa s > (x-d) (3457) 


where 
J — = J 
W. E(WIH9] 
/VAR 


d is the equivalent signal-to-noise ratio 


; E(WIH, ] 
VVAR 
M T. ; 
ema. Se (qt dt 
_ nel N 9 (54-48 } 
alg ae 
Jef a2(t) at +? s Q(t) fF au) Rg(t,u) duat 
0 0 0 
and o. (8) TS tne complementary erro! LUncCt 1 Onmewewcic: Kind 
. 1 co _y2/2 
o 8) —— fie dv (5.49) 
Yv27 g 


By oon a Gesdred Prova liltywort talsc alarm. x 
may be determined from equation (3.46). The probability 
of detection is then determined by solving for d and using 
emmation (3.47). 

PyesUbcrerterrac (5.4) Said (5.427) Sand usins (3225 )5 


fae numerator of d is 





7 =~ 24 | —— 
we YS : 
ee! 3 ee parang 


=, eu. Ff 





M T My | fe | 
2 Bie WE mGt ce edt = 2 ae e = = (Seo0e) 
meiy 2 it n=] a k a r 
2M k 
while the denominator of d (see Appendix D) is 
= 1, 
w fe 18,1 7? 
YVAR = =; | ft —~- (37505) 
Z| ae Nee 
2M k 
Therefore the effective signal-to-noise ratio is 
1M. ae fo tae 
d = M % a 2M Z — aM. , (G2 51e) 
n=] N k fea: ae i 


IPertedge of d along with (3.46) and (3.47) completely 
emtites the performance of the threshold detector. 
The equivalent signal-to-noise ratio for a target whose 


amplitude and position are known is easily shown to be 


1 


a at |: je? — I ae 
ni an satis KI 2M (3.52) 
n=] N k rc cs WN AL 


By letting the background be zero it will now be shown 
merce reduces to the equivalent signal-to-noise ratio for 
a known signal in white Gaussian noise which is well known 
[Ref. 4]. 

The energy dissipated during one observation in a resis- 
momeot | Ohm if the signal amiGte ox) is the voltage across 


miat resistor is 


oS 





T 09 CO 
Ms — > x, 
Beet, > fpr y ey) 2 ES (14) o5 (at 
Ga Kk UG 
2 > 2 
2a » Selves) : @62'5 3) 


time voltage resulting from M observations then is 


a ge = 
i nM 


M M 
” 
n=1 n=1 


oO 1. 
ee lhe : 
a E IE. (F,) | | , (aot) 


semester tuting (3.54) in (3.52) with en 0, 


d' = Ye f2M = /28 vm , 
N N 


which shows the familiar result that detectability depends 
Smeeupon the signal energy and the signal-to-noise ratio 


increases as the square root of the number of observations. 


36 


: 
; 7 a? 
oS se “ _ - 

2 Pe ee cle 

_ " 


>... e a {> 
(ere) ‘= : 
= ——— 


~ <— 


te 


7 Ss 
a 


, 
or 





— 
a 


Pee tec omOL De Nea rec TOR FOR NUTATING SYS LEM 


iimomeilierared detection system the output voltage to 
be processed is obtained from a detector-preamplifier com- 
bination. The irradiance on the detector Pi W/m 1s con- 
ferred tO vOlts with a linear scale factor. The voltage 


Sm@tput can be written as 


Gee a) = KH(t TQ) (aah 


where Elta) 1s the irradiance on the detector from a target 
eemosition cae and K is a scale factor in (V/W/em?) . LS O 


Miewecovariance function may be written as 
Rpt) = ELKB(t+r) KB(t)] = K°R, (x), (4.2) 


mimeres s(t) 1s the irradiance on the detector due to a partic- 
Meemmoackeround scene. The output voltage due to this sce: ; 
1S 


oe Re Gt) Co) 


Pt PpErOblcmmow 1s to relate the détector-preamplifier 
Output to the equations derived for the threshold detector. 

It was shown in Part III that an optimum processor for 
weak signals was a threshold detector that perforned a 
“emtain linear operation on the data. The linear filter was 
specified in terms of the eigenvalues and eigenfunctions of 
@ie intepral equation (3.18), where the kernal Ry (t ,u) is 


the covariance function for the background. 


By 


Biedrto Li the covariance function was noted to be 
Meriodic. Rewriting (4.2) in terms of the output voltage 


1 jkw, (t-u) 
Ra(t,u) = k E By e (4.4) 


Substituting (4.4) into integral equation (3.18), the solu- 


tion for the eigenvalues and eigenfunctions is trivial with 


Jkw ot 
_e 


b(t) 


and 


ee! 
dv a | BLT (op) 


where TJ is the nutation period. 
nee sicnal coefficients were also noted to be periodic 


and could be expanded in a Fourier series 


_ Jkut 
v(t) = K eee : (1e2G)) 


we Pas 


Recognizing the similarity between equations (4.6) and (3.41), 
the signal coefficients for a known signal which is expanded 


using the basis set above become 
2 Z > —. 
EMG = Aa lvt (4.7) 


For the threshold detector, the signal is averaged with 
: . x “ =e . 
Peewect to an a priori density function for Tq: Averaging 


Samepe done with the coefficients to give 
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7 (ype Ds (yale. 


kev of opr sid (4.8) 


Ky Tale 


Using equation (2.10), the coefficients ae can be written 





as 
| HRY GE ane axe ; es -jng ew 
eo far p(y) (eae) ee ye J, {20 fear? d“k (4.9) 


The transform of the radiance distribution from a target, 
Eek 57.) 1s a function of the target radiance, R(k), the 
mesrtion of the target, — andesite soot (Spreads sume troll Ob 
fmem optics, BOCK). Rewriting) (4.9)), 


270 \ Keke | a 


: ; Z 
A = 5 faekrce) t*C) J, 





-j2uk,r 


> > O ; 
far, p(t Je (4.10) 


femere density function, p(t.) , 1s assumed Gaussian with a 


Standard deviation oe then the averaged coefficients hecome 


I eee 
~2mp 9 (ky +k) Te oe 
‘ ; 270 4! k tie scl k. (4) 


x Y 





H, = 3" [R¢(k)e tA(kje IMO 


The Gaussian assumption on Ts is reasonable pecauee To 
often represents a pointing error from some designated tar- 
get position. Then Py describes the accuracy of the pointing 
iaemeanism. Substituting (4.8) into (3.42), the filter func- 


mr £Or the threshold detector becomes 


Vt 
? aan = 1 5 0 
ae) = he Hy v7 a ae (412) 
k Pa an /T 


So 





a oe re 





Rewriting the equivalent signal-to-noise ratio in terms 


Of equations (4.8) and (4.12) one has 


Mo : 
yo oF a, VT VM 
(2 a ep, Hii, |? — (4.13) 
YN7Z k i ek 
N72 


jmeevrlt be shown in Part V that By may be represented as 
5 cee 
By, = Op & Q, 


where ait is the amplitude of the background noise correla- 
iuomeeainction , 0? is the instantaneous field-of-view of the 
detector and Q, to Aico o ry ReLUnetleom OL “the Optical and 
HeteaetOr Parameters and correlation length of the background 


Meoraance . 


Define the following terms: 


1 M 
K M YY ail 
Signal-to-noise ratio (SNR) = Le (47 alearan) 
VN/2 
0 2 QKVT 
Background-to-noise ratio (BNR) = ——— (4. 4b) 
YN/2 


Deamee(4.i4a) and (4.14b), d may be rewritten as 


1 
@ 


1 


ee (4.15) 
1 + MQ, RNR® 


- 2 
d = SNR “WM |x |H,| 
k 


Noise e@uivalent irradiance (NEI) may be defined at the 


@-tectoOr-preamplifier output to be 


40 





ee 
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N 
+ + 2B = YNB, (4.16) 


Kee = 7 - 


where be i winemereetor- preamplitier bDamawidth . 
The SNR and BNR may now be written in terms of the NEI as 


NIT’ R aes 
/M SNR = hat = SNR! VNTB 


where NTB 1s number-time-bandwidth product and 


Oo 
/M BNR = ane /NTB = RNR' /NTB. 


Expressing d in terms of NEI one has 


wv 


[ 1 
2 


eet 52d Mae ee CL 
1 + Q, NTB BNR! 


d = SNR! NTR x [Hy | 
Ik 
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V. SOLUTION FOR RECTANGULAR DETFCTOR 


The equations in Part II for the detector-preamplifier 
Slept voltage and the correlation function of the background 
mus t be solved to obtain the filter function and specify the 
detection system performance. 

The two general equations obtained in Part II are 


juk-r 


ea = iaAe a On = ene 
H(t) = fuce B (k) t* (kJe 


Bei aie koek? 
n ony: = 


where H(t) is the mee Fourier series coefficient for the 


Padm@ati1on on-the spatial filter and 


D2 


2 2 2 2 
8, = firGol? [Fodol? wyGa Jo [2m fkgek’ ark 


~ 


where BL Ss the nth 


Coectfite@ent or che backeround correlation 
miiteeron at the spatial filter. 

A rectangular shaped detector was chosen for the calcu- 
lations because it exhibited the Simplest form of the 


saquaeions. For circular nutation the filter is represented 


auapmrecally in Figure 5.1. 
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rabcur ec Sseleme NitatlonewitheRectangudiar Detector. 


ferme detector function is given by 


1 x Wik een 2 
o/s ny enV 
ea yt) - = il . 
0 otherwise 
mene its transform is 
512m KX ot Ky Yo) | | 
T(k, sky) ms ee sin(k, w)sin(mk,h) 
ey 


Some ha t 


) 2 = Z 
Coee (wk w) Sin Cie) 
es Dee 
1 ky TI Ky, 


|x (k)|* = (5.1) 


k 


inempoOintuspreaa function of the optics 1s assumed to 


ies Oalssian: 
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|e 
7 2 
_ 1 20 
Rane i 
2710 
/mecscetranstorm 1s 
-202 6% (KEtKS) 


FY Ck) =e (orc) 


The correlation of the background radiance is assumed to 


Hemmer che form 


dp (x,y) = of e tlxl Bly | (525.4) 


where ae and Su are the correlation lengths in the x- and 
Meme ctions, respectively. The Wiener spectrum is the 


houmber transform of dp: 


9, 20 2B 


B 2 a (20k) 


=~ --; (5.3b) 
a + (27k, ) 8 


Walk) = o 
PeeitEicine (5.1), (5.2) and-(5.9b) Into (2.10)— the 


maatance from a point target becomes 


; a ae a 
Sy 27 Ciera tke ees ec den One araees he) 
BL (Gens = i” fe x x YY «4 Xe. 


J2m(kox +k y_) “sSanniow Ssinnk hi 
XO yo x Vy 


i) a ae 


-jn tan" (k,/K,) = 
e y J_lanp [ko+k* |\dk dk. (5.4) 
n a Bey ey 


Meeting 


S 
tI 


V2mTGK, 


= 
II 


V2m0k,, 
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equation (5.4) becomes 


a fie 7 (ugtus) 
GS) = ac 


TT 


2 
Pe [Cr -X9)u,t (ry -¥g)uy | 


-jn tan (u jae) /F 55 
C lea Mie | ales ae goat du_du Sat 
gee: x y 7 


>  , ; 

et as random, the radiance from a target may be 
averaged with respect to the density function p(T) as shown 
in Part IV. In particular, if a is Gaussian distributed, 


with standard deviation De the averaged coefficients become 


2 


eee : . 
re mf 2H Oo eS ee 
= 5 e e 
n 
sintk w sintk h =i tan” (k)/k,) 
Tk wk e 
yO eee ea Bo 
n ae A 2 
where 
Z 2 2 
a Fo 
oe en) 2 


tL : 
Crs V2mYK 


t= VInyk 
uy a, 


ooiation (5.6) becomes 


ey 





Wy ae ee = = eee 
_ Of (uy un) oy xXo*ByYo) : ji teat ee 
n L 


/20 2,22 
eh 4 pay euy Jenga (oe) 


Peco acinally substituting equations (5.1), (5.2) and 


mom anto equation (2.22), the correlation coefficients 


become 
a al Gael Stn ade J) 
: - fe Pt See 
Nh - > 
(wk, ) (mix) 
Zee Z 9 
of ae ey seep ae Jn 2710 [ke +k° dk dk, 
a® + (2k,)° BY + (2mk,) y 
(5.8) 
Nor ting 
ul - 279k. 
uy = 270k, 
eomation (5.8) becomes 
9 
cy en 2 era uy? sin = ut : 
Meteo es | ey | 420 | 420 | 
B. = ——x~— |e = 
n Z W tt tt 
Yee ey ee a ad 
20 /20 > 
Tie Wa io) \ 
(ag) “tuy (Bo Y+uy Y 
where 
2 
Q¢ = wo : (S29) 
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Mer dntecrals ian (525), (5.7) and (5.9) were difficult 
to evaluate. Although the components of the integrand are 
relatively re functions, the argument of the Bessel 
m@metion prohibited separating the integral into two one- 
dimensional integrals. No closed form solutions were found 
so the integrals were evaluated numerically. The numerical 
method used is based on work by Pierce [9] who applied 
Gaussian quadrature formulas to two dimensional integration 
by integrating over a planar annulus in the (x,y) plane. 

Gaussian quadrature formulas are a means of evaluating 


mmemantesral by summing weighted values of the integrand at 


ieeiii1c points. For the one dimensional case, 
b b 
orc anaes Wx) ft (x )edx sarki} 
a a 


where g(x) is the integrand to be integrated and w(x) 15 a 
weight function for which the snecific integration was 
derived. . For well behaved functions the integral may be 


evaluated as 


b , 
i wee (x)idx = 
A i 


lms 
bem 


A f(x; ) ser Om (Sait) 


where the X'S ane the mezeros Of Ene mee polynomial 


P(x) = WM (x-x,) 


loos 
4 


of the set of polynomials mutually orthogonal on the inter- 
moe [a,b] with respect to the weight function w(x). The 


weights are given by 
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Ne =F w(x} L; (x) dx (oom 5) 


where L. (x) = 7 COPAG Cer Sp Lede coe is the Lagrange inter- 
polation coefficient. 

For a weight function w(x) and interval [a,b] the poly- 
nomial POX) amcderts Zeros Xs and weight factors A; need be 
computed only once. For a number of weight functions and 
intervals the set of orthogonal polynomials is known. Stroud 
and Secrest [13] give x. 'S and A.'s fOr sae Varvets. On elon e 
functions and internals. The degree of the formula deter- 
mines the number of X= 'S and A.'S. A 2M-1 degree formula 
Piienmave M points and is exact for polynomial integrands ae 
@emnree 2M-] or less. 

Pierce applied Gaussian quadrature integration to two 


faoaemsions where the solution is of the form 


Nei 


Cd 


ffe(x,y)dxdy = =f Doers ec 
1 J 


The integration is over an annulus in the x,y plane wi h 
inner radius R and outer radius 1. Rewriting the integral 


in polar coordinates 


wa: 


i pC) Cremels Miceli ciel ssue-aey ae) ae Derr Gh Geers) 3, | (Sc) 
0 oe) ae 


f 
R 1 


Pierce showed the summation could be rewritten as 


4(m+]) m+l - 
y ee Gi eC OS oat sees WC ae 7 Beet 
1=1] j=l aa * it J 
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where A. and Be dre Cie welts ton wre radial and -anvular 
mirections respectively and 4M+3 is the degree of accuracy. 
For an arbitrary annulus the formula may be obtained by 


femrbtang equation (5.14) 


Zia <2 
I = f f ie ee recaloctal c) (5.5 1G) 
0 Yr 
5) 
pod betting 
= 
Pp a 
3 Ze 
I = r5 f if of (r,p,9) dodé. (Sele) 
O a 
l/r, 


The approximate solution can be derived as 


Cant ly mt i 5 ( 7) 
T= ¢£ UAE pina 9 Sree (cs CMe Se ie ge eons} 
<a j=1 iu) ewe) \ (oe i 
where 
A; = 2u/4 (m+1) 
and 
B. = (Geese ye w. = weight 
J et J 


Memeers the order of the orthogonal polynomial, in this case 
exemLegendre polynomial on the interval (0,1) and; ''s are 
@eemezeros of this polynomial. This type of formula is known 
as a spherical product Gauss formula [12]. 

The annuli were picked at radii corresponding to the 


Zeros of J, (2) because J, 2) is the fastest damping component 
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Sine Integrands. One additional annulus inside the first 


memonwas found necessary to account for the term when 





atx 
@ 1S small and x (and y) approach zero. The distances 


between annuli outside the first five to eight zeros were 
mumicenot to be critical and thus spaced arbitrarily. 

The Gaussian Quadrature formula used was a Gauss-Legendre 
24-point formula as listed in [13] thus giving a degree of 
meemracy Of 99. The integration 1s limited to only the first 
miecrant aS mentioned below. Thus for a 24-point formula, 
(M+1) 4 or 576 points per annulus were used to evaluate the 
immeeoral. 

Mie expression for i VsSeecashiyeseci to WC aAnveven mune 
Prenmmoocth radially and about both planer axes, therefore, 
only integration over the first quadrant was necessary. 
Vemma little more difficulty, H(t) and H can be shown 
(cer ppendix FE) to exhibit the same property provided one 
expression is used for even n and another for odd n. By 
Calculating only over the first quadrant computer time was 
considerably shortened. 

Inmmorder to solve for coefficients of arbitrary order, 
micempesse] functions of that order must first be obtained. 
The available library subroutines were Rommel ines Sewince for 
this use for two reasons. First, these routines have an 
memer limit of 100 on the order of the Bessel function that 
eae be Computed and second, the routine had to be called for 


eer individual order. 


=) 





ime recursive relation 


on 


GY = Fo) - Iya) (5.19) 


Mmievrdes a basis for generating an array of Bessel functions 
of various Cudeverm  HOWweveT tor accuracy J+) ) and ae) 
must be known. 

The most accurate method to numerically calculate Bessel 
functions of various orders and arguments was found to be a 


Maeeorm asymptotic expansion involving Airy functions [1]}: 








4 273 
; A.(n?/Say © ay (A) 
J_ (nz) = siti : ” 
n 1-2z as ie 2K 
n k=O on 
at(n?/Say by (A) 
emigre 0) 
n k=0 on 
where ae and A; are two types of Airy functions and 
= asl? = in| 1B Vie? |/5| - V1-2¢ 6 Il 
a 4) = ] 
ee) = 748 Ae ae ry + ae | 
Ae OC (Oey a ene to eae 
a, (A) -sS a (A) 
camel (S220) 


by (A) << by (A) 


Pie Airy -tunctions are calculated using the equations 


ail 
4 


1 -6. 
A. (x) 5 xX e- i {-36) 


and 


ae 
Ai (x) > aaa « eae 


ti! 
J 
NO] — 


a1 








where 


Bi 


O 
1 
WI 
va 


(Se22) 


anid 


x = noe A 


ime functions £({-6) and g(-6) are tabulated and linear 
himempOolation is used to determine their values. <A is guar- 
miimecad £O always be positive if n is selected to be larger 
than the argument. 

By picking an order much higher than the argument, v, 
and Ja using the asymptotic expansion, the recursive 
mourceronship (5.19) may be used to generate an array of 
Successive orders of the Bessel function down to and includ- 
ing Jo: However, error builds up rapidly using this method. 


pee cond relation 
1 = Jo) t 25,0) +r 2540) Ess. ace (S272 5)) 
Maye Used to generate a normalizing factor 


C= 1/(Jo(y) + 23,7) + 24) + boa | (8) 


Miikerplying each Bessel function value generated in the 
recursive Pea anien DY C resulted ai very -aecurate values 
being obtained. The computer program used to solve the 
fam@ations is listed as Appendix H. 

To evaluate the numerical results a check was used which 
summed the coefficients and compared this summation to an 
analytic expression for a sum that could easily be solved 


numerically. 


a2 


a = eee ee 





meee pareicular phase angle of the nutation cycle, Wot, 
trem adiance function for a target is 


Jkw ot 


H(w t) = Here (5.25) 


var lees 


mea point Carget and a rectangular, circularly nutat- 


ieee CeCtor 


-2no°k* sinnwk Jémk, (pcosw t - (r,-x,)) 
H(w.t) = fe * —,~— e i 
O es. 
a), a: 
vireo Ko Slik 1 2k (pS = (rey yy) 
ae y ee aie 2 OS cl 

yy 

(5226) 


which has the solution (see Appendix F) 
idl : ( : 
H(w t) 7 Fy (ry x90) Py Cry x95 || Fy Ory ay gh) 


a E. (r,, Vo i) 


where 


r-z-pcosw t , 
EB, (v,2,0) = | <a 
/20 2720 








r-z-psinw .t ” 
E,(1,z,a) od eae + j 


V20 2720 
pier t 1S the error function 


x Z 


Erf(x) = fe en ee ec (5.27) 
| 0 


|» 


TT 


Byeevalvating thesanalytic expression at. the phase an 


stant (wt) that the detector crosses the target and comparing 





Memewainst equation (5.27) for an appreciably high n, a 
reasonable check was made on the accuracy of the individual 
eaeomal coeificients. 

Mer the coefficients averaged over a 


, the check was 
made at wit = 0 where 
CO 


H(o) = E Hy (5.28) 


aiemecie analytic expression is 
Etro [Ey (ryoko sw) - By (rx -w) |[F3 (ry sy gh) 
4 rex 2077 a cae ecko “3 y?o 
Es (1, .¥5,h)] 


where 


E,(1,2,0) = Fo (1,250) (5229) 


Ww t=0 
O 
The rms background voltage was used to check the coef- 


ficients in the correlation function. The mean square vol - 


age 1s 


ae E{B-(t)} = R,(0) = (5.30) 


p ( 


ast S: 
WD 


Appendix B shows that the correlation function may be written 


as 


0 es 2 


3 Cae : , O 
Rp (Tt) = S| (k) | Wp (k) a tJ SN d 


5 k (16B) 


pies tuting (5.1), (5.2); and (5.36) am (1608) and evaluat- 


ing at t = 0 
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aes 
9 sin’ (wk Ww) i : Ky 
Cp = Op2a {[—_———*, ay aS dk. 
(Tk ) Ca (27k) 
aed. 
sin’ (mk h) gant a0 Ky 
‘0,28 f-———~— _ +, dk, . 5505 


2 2 
k y 
(1 y? aes ae 


It can be shown (see Appendix F) that an analytic expres- 





Son 1S 

Bor nce R ( wiamnin 2) 

B RB b] > 
where 
(a 
R(u,a) = “5 )1- Erte w) 28 has Be 
O 
et YT 











LE 
oo 0 es ‘ 
a e ie Erfc |jao + =< 
3 (iio) 
20. 
+e 4 Erfe fao - = = 2 rfc (a0) S52) 








aideerrc 15 the complimentary error function 


00 2 
Erfc(x) = ee ed (5.33) 
x 


<a 
v1 
The computer program used to solve equations (5.6), 


fone) and (5.9) and calculate the infinite summations 1s 


listed as Appendix H. 


Sys) 





Vie NUMER TCH ial is 


homerical results constitute solving for the signal and 
background coefficients and using these coefficients to 
specify the form of the filter to be used in the threshold 
femeetor as well as to calculate the probability of detection 
feomaetunction of signal-to-noise ratio (SNR), background-to- 
noise ratio (BNR) and number of nutations (M) used to make 
M@iemdecision. The probability of detection is determined 
for the averaged equivalent signal-to-noise ratio and for 
oecttic point targets. 

Many different values for system parameters and target 
iieeerons were used to compute the coefficients. A typical 
Setunot Parameters ist; 

itieaeron radius — 15 


width of detector = 30 
Memon ot detector, — | 


blur vemrelesstandard deviation gam. 

background correlation length x-direction = 20 

baekGromid correlation Jengrn y-digeetion = 20 

MOSTLiOn Of detector im detector eeeomdinate sysecem 
CEdsreet ten. =) b> 


Vesciintee @ te lOllenn an 
target coordinates (if desired) 


dire creron — 
y direction = 0 
POoimecIineg Error Standard devia tine an, 


Peeagletector utilizing the parameters above would trace out 


the area in the image plane shown in Figure 6.1. 


Iaq nt t>oocrve -inema Lapadians- 
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44st 7 7+ 


re eee AY 
(WML LIL Lhe 

t = 31/4 

Psineo. iW scanned Area OL @iitat in etcetan. 


Tmis type scan was chosen because it crossed the center 
Gameme coordinate system which would be the center of the 
meeret blur circle for zero pointing error. Also, in the 
meteor 2 random pointing error in thé absence of bias ze > 
teers tne center of the coordinates. This 1S important 
because the threshold detector design is based on a Gaussian 
meeting error with zero mean. 

Aesample calculation of 121 coefficients using the para- 
meters above is listed in Appendix G. The equations shown 
jorPart V were used to check the accuracy of the coefficients. 
iipbe G21 lists typical values for the summation of 121 coer- 
pierents fOr various sets of parameters and the calculated 


Summatzon values. 
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TAB EE Ont 


SS 
oe — ee = min : - : = 


Calculated Sum Actual Sum 
Saemal (point target) .6826894 .6826993 
.6826895 .6826994 
-5438002 -5438087 
Signal (averaged) .2990664D-6 a7 Oe 52) = 6 
Cole se D4 2 525g 21) 8 
eZ oles Jom aa 20 O15 a 
Background .6013477 woO 72.50 
~2347420 $2 5474.20 
NOUS .8319764 


Besides the general shapes of the background and signal 
miectma, Variations in the spectra due to detector size, 
Meaekeround correlation lengths and coordinates of the detec- 
mMommare of interest. Curves in the following figures show 
the effects of these variations. 

Memeuires 6.2 through 6.5 show the envelope of the one 


side of the two-sided background power spectral density 
Sp (w) = : 81.6 (w- kw) 


where WA iis Mormalized towlie s (hee cUnVes (Ol Pivures O47 cmd 
memetinstrate the effect of varyine the background correla- 
meagmetencths. In Figure 6.2 the correlation lengths are the 
eameein both x- and y-directions. Practical measurements, 

however, indicate that correlation lengths in the horizontal 
fe) direction may be longer than in the vertical (y-) diree- 


mieme ine curves in Figure 6.5 are for this case. 


ons) 








ie scurves ot Pigure 6.4°and 6.5 relate the spectrum to 
@aemeees if Size of the detector. In Figure 6.4, only the 
detector with BNR equal to 9 and number of nutations equal 
momen are Shown in Figures 6.11 through 6.13. Each figure 
shows the envelope of one side of a double sided spectrum. 
Two cases are considered. First, threshold detection systems 
designed for different Gaussian pointing errors but with the 
Same Size detector are shown in Figures 6.11 and 6.12 and 
second, detection systems designed for the same printing 
Seuvor DUt utilizing different size detectors are shown in 
Pemomes 6.12 and 6.13. 

ie curves Shown in Figures 6.14 through 6.16 show*™the 
Wagemed filter to a particular point target located within 
tiemsean. Notice each of the filters exhibits a band pass 
@memeeteristic. This 1s because of the additional noise 
eemmenent (background). The oscillations observed ain the 
MOint target spectra of Pace: 6.09 tO 6355 have beenwcup- 
messed to provide a look at realizable filters. 

mares probability of detection for the threshold detector 


was shown in Part IV to be 


foe) , , ] L 
Q, = SNR YM} & el aie en) 
; 1 + MQ, BNR 


oe 


raeures 6.17 through 6.20 show QG 25oa Tome t lon Of SNR, 
BNR and M. 

line probabilities of detection for point targets located 
a@uespecitic points in the scan have also been plotted. The 


probability of detection for a specific point target is shown as-the 


ste 





—— - > 





fae OL the detector is varied. AS width is increased, so 
/mmene area scanned by the detector. Figure 6.5 shows the 
mee rence in spectra between a, rectangular detector and a 
Square detector which have approximatcly the same surface 
area. It should be obvious however that the rectangular 
MemecuOr SCans much more area per nutation than a square 
detector of identical surface area. 

The envelopes of signal coefficients (magnitude only) 
fogera pOInt target are plotted in Figures 6.6 through 6.8. 
Phase information depends only on the location of the target 
relative to the initial point for the nutation cycle which 
is along the x axis. This has not been plotted. Changes in 
the magnitude due to changes in detector size are shown. 

Eetnes 6.9 and 6.10 show the envelopes o£ coefficients 
for the signal that is averaged with respect to a Gaussian 
pommm~ans error. The bandwidth of the averaged signal 1s 
much less than a point target as expected. The figures she © 
the spectra for two different standard deviations of point- 
ii error. 


Matched filters (magnitude only) for the threshold 


Poe 1 
Sep are ee Ea 
tr k 1+MQ],BNR2 


Q} = SNR'VM . 
= 2) 1 Pa 
rlé|? —1_, 
k 1+MQ, BNR 
where 
1 M 
SNR ~ oy 7 a vt LS 
n=l Vi 2 


60 





~~ 


y™ 





Figures 6/21 through 6/24 show QA Ase tive @erongo br ao) ae 
BNR and M. 
fare parameters of the nutating optical system have been 


mepeed to show their effect on Qy and Qi: 
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VII. CONCLUSIONS 


The problem investigated in this thesis was to determine 
an optimum processor to be used in detecting targets with an 
infrared nutating system. The nature of the system divided 
ine problem into two areas: optimizing a spatial processor 
according to size and shape and optimizing a temporal pro- 
€essor that takes the output time varying voltage and decides 
if a target is actually present. 

memeral equations for the output from a nmWtating detector 
were known from Samuelsson's work, however these equations 
had not been solved for a specific detector. At the same 
time it was observed that the form of the covariance function 
moeetie noise offered an easy solution to integral equation 
of the Karhunen-Loeve expansion which made statistical 
detection theorv appealing. Some work had been done on tr: 
mmm Of temporal processor uSing Statistical detection the sry 
but this was limited to only one nutation because of the 
common background noise component between nutations. 

Harger's derivations provided a means for describing a 
memperal processor using statistical detection theory that 
besed its decisions on multiple observations and included 
eme background noise. His work was extended here to include 
unknown parameters, amplitude and position, which more 
eiloscly characterized the system. This extension, called a 


mereshold detector, was shown to bev optimum in the case ot 
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morals signal-to-noise ratios where detection is most diffi- 
@uit-. the derivation also used the Gaussian assumption to 
@escribe the background. This may or may not be correct. 
Poeein practice, the detector's performance may not suffer 
greatly if this assumption is wrong. 

To implement the threshold detector, the integrals 
describing the detecor output were solved using the Gaussian 
Quadrature method of numerical integration. Checking the 
computed coefficients by a summation provided a means of 
Weeepcing the validity of the integrations. The form of the 
mopeshold detection system using a rectangular detector was 
feeermined and the frequency spectrum of the optimum filter 
Was shown. To specify performance, the probability of 
detection is plotted against signal-to-noise and background- 
tO=moise ratios and the number of nutations on which the 
decision was based. The signal spectrum for a point target 
and background power spectral density is also plotted. 

A designer may find these results useful in developing a 
System or investigating the performance of an actual system as 
compared to the optimum. He may also use these results to 
@eeermine the size of a rectangular spatial filter to be used: 

Ricwre work should ancliude extending (these Tesimits ire 
aeeial filters that are circular or elliptic. dhe same 
melterical integration algorithm could probably be used. The 
threshold detector was found to be highly sensitive to back- 
ground correlation lengths. Thus, some form of adaptive 


meocessor should be included to minimize the exrfect of 
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background. At the present time, a measurement program is 
underway at Naval Weapons Center, China Lake, California 
memaetermine average background correlation lengths but for 
different environments, these correlation lengths could be 


expected to vary considerably. 
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APPEND Lee 


DERIVALILON OF SIGNAL SPECTRAL. COl er hernia > 


live radiance on the detector 15 periodic because of) tie 


Nutation therefore it can be expressed in a Fourier series, 


oo Jnwt 
oC) resi ls ane: (1A) 
Nn 
where 
A = Vitae 
N O 

H = 2 | H(te dt (2A) 

vine 
; 2 
H(t) = JN'(r)e(r-p(t)) a’r . (3A) 


Em@iestituting (3A) in (2A) 
= Wise 


T 
See ' ue : O 
ne = fe r N (3 7 J dt ple o(t) e (4A) 


The transmittance, t, may be expressed in terms of the 


inverse Fourier transform, assuming an infinite image plane, 


2 -j2nk+ (r-0 (t)) 
r(r-p(t)) = fa’k tke . (SA) 
Pevwriting the exponent, 
oie) = k-r - pk coswot - pk sinw,t (6A) 


Equation (4A) can be rewritten, 
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. T -j2 ker 
+ fy ] Diao po aie 
nay eet eOs) ri dt fa kt*(kje 
jy2mAsin(w t+6) -jnw_t 
a O 2 O 
where 
A = ke + ke 
x y. 
6 = tan + evel (7A) 
xi “y* 


Observe that 


9 = 2 
N'(k) = fo rN'i(rje’ 77 (8A) 
thus 
T = J ioe teal) 
H = fa’ Nick) T*(KVeIME LL at i. 
n _ : ia 
JnAsin(w t+6) 
“e (9A) 
Weame the relation 
1 i VO! 52 25 
J_(z) = | ee dee do, (10A) 
n T ie 
~memexpression for the coefficients becomes 
7 2 ; ok jné ie 2 
H_ = fe k N'(k)t*(k) e J, [2x0 kotke |. (11A) 
et 
‘ ge ee 
@ = tan ky/kx 
then 


jee 2 ene 


Rewriting (111A) 


&9 





_— -_ a ec 5of wor. 





eal 2 : ’ -jno Z 
HS j fa k N (k)t (k) e J, [2m [+k (12A) 
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APPENDIX B 


Dery lON OF BACKGROUND CORREEATION COREE) CEI ils 


th 


aire 7 power spectrum component of the background at 


fee auency NW, is given by 


ampeny 


: O dt (1B) 


T 
a = ee Rp (te 


But Rp (7) 1s the correlation function Of ede tecetor soulmouE 


mieraced with respect to the starting time. 


Rp (1) = 7h E[B(t+r)B(t)] dt (22) 


Beame the expression for B(t) from Part Il, 
leat 2 
= 1 = 
Ra(t) = rl dt E[ Jn COGS (eels 


2 
JN" Cry) <(r-8(t9)| 4 a oe 
Interchanging the order of expectation and integration, 
R -1f' at fa*r, fa’ r,-e(t+t) | 
Bee ea Te aati Es aio ads 


-t{r-9(t)]E[N' (ry )N' (29) | (4B) 
Assuming stationary background, 


E[N'(ry)N" (12) ] = @ (E722) ($B) 
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ee 





where Op (41 7Ta) iS the backerommed covariance fUncEeton “Onetiic 
image plane. 

ihmissnow convenient to express the intesrand ins terms 
OG two-dimensional spatial Fourier transform. To begin with 


WALD Seat 7 


t(r) = fr(kje ~~ d°k (6B) 


where t (k) is the Fourier transform of the aperture function 
TAGS) and k = (K+ kL). Then 
15 [r4-¢ (t+t)] s [r.-p(t)| 
= fa°k yr (ky) fa2x,ct (kyon nd T2272? 
-j2n[ky 9 (t+t)-k0(0)] 
e (7B) 


Smipotituting (5B) and (7B) in (4B) and interchanging the 


Geacer of integration, 
Z 2 
Rp(t) = fa kK, ky) fark, t*(k,) 


f i -j2n[k,-p(t+t) -k,*9(t)] 
e 


OT: 


mfr 


j2uk =) ahs 


. soe) & Bn 
Nice Teel { 2 wea? Ko 25 
: fa Tr, © dr, op(t,-rale 


(8B) 
Iiteslast integral becomes, with R= T17To> 
2h Wale) COUN ed 
| fa?Rer re fs “| iw (9B) 


But the bracket above is the Fourier transform of a(t), or 


the Wiener spectrum of background on the image plane. Let 


WOU 2 
W2(k,) = fog (Re a“R (108) 


92 


— er U7 — ss. re ae oe a 
yA 7. 42 2 





If cs} VS the fourier Lranstorm on the porntsspmeda wine 
tmen and Wp Ck) 1s the Wiener spectrum of background on the 


momect Plane, then 
trey) = 2 
Wak) = [FGCK)I" Wyk) (118) 


Mims relation is analogous to the output noise spectrum 
Eyeeression due to a noise input to a linear system. In terms 
ONE We (k) tebe Ae tiem Oto bass 


g 


erie 
K) fe ry 


We d = Wi (k)6(k,-k,) (12B) 


Semesctitutine this into (8B) and performing the integration 


Weegee respect to k,, 
2 Capen 
Ry(t) = fark [x(k 1? WE) 


D 32a (oUt) eee 
i e “oS ~ dt (13B) 


re 


O 
Eemee p(t) = (p coswot, psinw t), 
k ~ = + + 1 + 
k [o (t+T) p(t) ] k,pcosw, (t T) kee ina. e)) 
“k pcoswot - Koes” 


= A(t)cosw t ‘ B(t)sinw t = c(t)cos(w t-v(t)) 


where 
A(t) = p[k cosw ot + Spray 5 kV] 
B(t) = e{-k,sinwt + kK, COsw at - Ky] 
c(t) = fA? (x)+B7(x) = p [2 - coswgt) (kg+ky) 








and 
a dee ey) 
p(t) tan A(t) (14B) 
iimemetore the last integral of (13B) is 


12 J2mc cos(w tty) 
rJ e dt = J, (21) GES 


where Ji0) 1s the zeroth order Bessel function. Therefore, 


mmemcorrelation function Ry (7) 1S Piven Dy 


Ra(t) = [| (k) |? W2(k) J (4npk sin —2 lack (16B) 
where 

ee eee ard) ee eee 

an  & y re ~ xX ~V 


Pinally, the power spectrum coefficient Bt 1S, £romel le) 
= ee tf Z 
8, = Jit’ wE(Ka"k 
a WT Sq eu 
; | 5) (20 Sin 2 | e sks (17B) 
T : O 2 
where a = 2™mkp. The last integral becomes with x = Ww ot/2 : 
Lie 
ah el Can sint)eCOos 2m sax 
ey O 
1 f 


- j= Jo(2a sin <) Sine wddx (18B) 


Mire cosine integral equals 3“ (a) and the sine integral 


vanishes [2]. Therefore on using (18B) 
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B= JITG)1? [FSG 1? Wek a2 (2mpk) a7 (198) 


This is the general power spectrum expression for a nutating 
@emectOr output. The power at frequency NW, iS expressed. in 
feos Of aperture transform t(k), Optical tramster tunecion 
FAK); Wiener spectrum of background Wp Ck) and nutation 

radius through ceo The integration, 15 90Venrrumewen time 


Feooilane . 


eis) 








APPEND Ee 
BAPE CT Aron Or My 5 


ire essential calculation required to-find 


Ay 2({Z,3) = EptA,,({Z,}}B)1 ae} 
for Gaussian background B is the expectation 
ek at 5 
I = Ey jexp , ACtIBCIdt = cy Cea iteiae 
O O 
where 
c = M/N 
and 
5 M 
A(t) = y OE y(t) - Sy(t)) (2) 


B(t) may be expanded in its Karhunen-Loeve (K-L) repre- 


sentation provided the mean square value 


E(B°(t)) < @ (3C) 
and thus the integral equation 
se 
OV arbe ar (4C) 


AgOg(t) =F Ry (tou); Cu)du 


can be solved. ds are called the eigenvalues of the equation 


pals. (t) are called the e€1gentunctions. 
i 


the veorre lation. funect reno enc 


Using Mercer's theorem, 


baekeround may be expressed as 
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Rp (t,u) = Ay, (t) o* Cu) oe) 


ite K-L expansion of B(t) is 


co 


B(t) = 2 by oy (t) 


i 


with 
dt 
bh = SJ B(t) of(t)dt UE Sol (6C) 
Gye. 


where the coefficients, bys are uncorrelated and because 
Meee 1S Gaussian the coefficients are statistically inde- 
pendent with means Uy and variances dy. 


Also expanding A(t), 


co 


pac jl a, 1 (t) 
k 
Waecn 
4p 


ay = A(t) $#(t)dt, Ve (7C) 


Rewriting in terms of the expansions 


(i co 0 co © 
1 
i= 5 exp! V{= 2 2 a DE = 5 2h ae bl) oe Cee 
B | F ie ee k eee k k n 
T0000 f 
* 
= 1c Z a b bE o, (t) on () dt . CSiG) 
Using 
T s 
Sxy = j b (toy (t)dt 


of 











- Ne ee ee 
I Bread jex E (5 abe + 5 a*b_) a : ob (9C) 


Using the independence of the coefficients 














Ea ie epee ee 2 
I I En @XxD E aa + 5 anbe E ba | 
n 
aan ais a 2 
“wie xD ee exXxpaG a a @e es 
n 4c 
Por Gaussian random variables 
exp| w m,m&/(1-2wA)| 
BROS) 2B eas <a 
where 
te Eee). 
= = Kom 
r Eee m,) (x my) 1 Cie) 
and w is a constant. 
Let Cape cy 2c 
n n 
and ae is as 
then 
2 
aan exp {-clu,-a,/2c| (1+2er9} 
iE oxen One 7. = eee (2) 
bn n Cc Ge roa 








Smostitutines12C) anto.(10€), and sevagaeimestems 
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[2 Pe 
n n 
n 2 


: “1 1° 
—— II lh 7 2ch.) Exp = T+ 2ch_ 


ro) % = 
exp 4-= & Y Bn {Mn @n/o) 
2 A: 2cr, 
c N ry ipl 
a |. Pa eS So 


This result can be written in closed form by using the 


@eainition AL = (a,,) 
2 
ea Ba) OT 43 | 
x 2 4 ee 7M f dt,A(t,) i dt,A(t,)q(t,,t,) (14C) 
n n O O 

where 

GG, U a = 2M 5 ibs. b(t.) oe (t,) (15C) 

be eZ N 5 1+ 2c, 9 names Rise Vinay” 


meesonen that operating on q(t,,t,) with 


Teele 
i es Ry (t,,t,) + 
‘ome(e 





Ry (Ty 5 t>) (-)dt,dt, (16C) 


and using Mercer's theorem yields 


i 
a(tysto) +L Ry(T pots) a(tssty)dts = Ry(Ty.t,) (170) 


P| 2 


M 


Likewise, 


* = at koqgk 
ee ea) We Cuk-ak/c) 


lt SPA Le) 
n 


int jot (18C) 





ft 
=~ 1] fa tad Att} 
ay J my, (t) = 


@) 


3), 











where 


_ n 
A(t) = 2 -yyez7 6, (1) GEIS) 
n n 
such that 
ib 
Nh(t) + 2M j Rg(tz,t) h(tz)dt, = Mm, (t) Gua 


iaring the logarithm of both sides of equation (13€) and 


using the equations above 


an I= 5 (1 + 2cd,)7* 


n ,! 2 } 
Zs j dt, v2 CE ne 54( 9] 


— 


T F M | 
f dt. 4—. 2 . (2. (t.) = Soe cea 
a 2|N wap mh m* 2 ie 
T n 12 M 
oe CO ee antag 7 (Z(t) — SaCeiee leat hide 
. -1 
ene | ee Oe: 2bi) 
Tl 
eh oe M M M 
ee an oe ee ene Ca) Mem OS Gea) Oe ic) 
Ne a eee Coe eee ce le eae 


M M M M 
- ¥ S(t.) E Z(t,)+ 2 S,(t,) ES (t,)} q(t,,t,)at,dt 
Ae ne ey Be ee ee ae 1 aes alge 


- 
m, (t)h(t)dt 


' 
or. 


M 
Z (Z,(t) - Sp(t)) h(tjdt . CZ) 
ail 


ri 
\) ee 


Jnr 
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_ >» | — 5 = 7 *% = * ‘s_ = 





times likelihood ratio desired is 


9 M T 
A,2UZ,}) = exp jy mal : 2 (t)S_(t)dt 
i, ne 
mes Sa Cust) aie (22C) 
n= Ec 


Taking the logarithm of both sides and substituting the 


even ression for &_ CD; 


2 M sh 
Mei (zt) -5 2 sf 2 (tS (t)de 
n=l] o 
M T 
-iL x s s@(tjdt 
N n=l 0 
~ 2M | 
$7 1] + —A 
: inti 7 i) 
1 T T M M 
aay ae | Nees ee ee) 
NM — a | oe es 72 
M M M M 
= ¥ Sy (tz) =]. Za i ae ee ee eee 
g=] ae | a4 m2 ee WeeeZ = a 
M M 
+ E Solty) 2 S(t) Ce a) clay ate, 
Q=] m=] 


ik 5 T 4h 
- f m (t)h(t)dt on f i (Z, (t) - S,(t))h(t)dt 
O O =] 


— 


(25) 
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APPENDIX D 


KARHUNEN-LOEVE REPRESENTATION OF EQUIVALENT 
SOIGNAL-TO-NOTSE RATIO 


The variance of the threshold statistic 
M f - 
We= Z(t)q(t)dt 


where 


7 (late 
a(t) = 2 fice) - J Fapatu,t)au (1D) 


O 


under hypotheses H and Hy 1s 
nm ft -2 2 : - 
VAR = = | Get ate a fat-act) Jawa (uy Ry (tw) (2D) 
O 


substituting (1D) into the first term of the variance 


eepression 


it poe : : es 
ae q‘(t)dt = 2M af Pee) 280 | f(ujq(t,u)du 
O | O O 


eit : 
+l f f(u)f(v)q(t,u)q(t,v)dudv (3D) 
oe) 
Using the Karntnen-Leeve representation 


on COCO. T 
NM = 2 es 2 M * * 
3 q°(t)dt = = F Ce, J o, (t) o*(t)dt 


CG Coin a Bp. N - Als Ty 
tee) fectceeeee seen eee * j * 
220 REE oJ a (eveacevat fo, (woz (uddu 
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2m 2m 








Prnn fre Nm N “n e  (t)o*(t)d 
+ —_—_— SSS eee a 5 ig 
k mn Krk y+ on dr p+ 2m LAG Ng ue 

N m N n 


ie sh 
Jo ecw eg cu)ar J do (v)ozQdt] (4D) 


Using 
f 
= rs 
Soe J oy (Hegceyar, 
then 7 2M : 2M ; . 
NM -2 2M me Sats N k Nak 
— q (t)dt = = 2 f,f% 11 
“a Ny kk  . 1+ 40, 
ee ae | eee 
= FEE 7 aM, (SD) 
N ok 


The second term may be expressed in a Similar manner to 


give 
ae 2 De. 2 

van = 2M 5 pee f__ 2", (2M\" 5 bee 1 | rk 

Me Sse as ee Sle 28, 

N ok N “k/ 

We = 2 Be eee (t D) 

Se 

| N *k 


fhe Cauivalent signal to noise ratlogiine- ame Ils 


Me 2) Ales.» 
raf F(t)a(tat | 
ee Ape OR Ok we A rams (7D) 
¥ VAR 


m@ne K-L representation of the numerator 1s 


tome 
INO 


iy M CO ‘ 
2 i . 2 2M 
zi | se (rote Gc) alias ba y | £, | (i+ s,)] (8D) 


n=1 O k = 





o— Se he eee ee 6 —a = 2 


d may now be represented by dividing (8D) by the Square root 


of (6D) 
M - 7 2M 
Ze Sahl tl LLP ee 
n=] 
oa I 
1) ree, 2M 2 
ELE L 7 2 tap M) 
k 
eee - 19 mM. |? 
ME, an MELE / Ota Ay) 
G-= et LK... ae (9D) 


YN/2 





Pe ae ew PO AR OU ys a 


APE EN Oe 


NON-ZERO TERMS FOR NUMERICAL INTEGRATION 


meuation (5.5) tor the radiance function to1resaypome 


target may be written as 


I fe fox -¥ 9 9 EPx*Qy) Sinax singy 
n es Py 


tan’ 1,(¢ x+y? ) axay 


where 
a = w/v¥2o 
8B = h/v2o 
G= V207G 
_ 2 
|e = (Gea) 
_ 2 
Ge= eae 
and | 
> = tan! y/x (ley) 


WaeE(x,y) be that part of the integral that 1s even about 


both planar axis which is 


0 2 2 . . 
J oy Oat rey 
Oe) ts 7 © mie = s3ey J (« Xx a ), (25) 


TT 


Using Euler's relation the integral becomes 


EOS 





1 = [E(x,y) [cos (px+Qy) + 7S5in( Px+Oy ) ) teosmd. — 15 iam Iced, (SE) 


Using trigonometric identities 


= JECuy) x TT. dxdy 


i=] ~ 
where 
eT = ~COs PX cos Oy Costing 
T, == Sine x Sim Oy Cos ta 
T. = -sin Px cos Qy sin no 
Ty = -cos Px sin Qy sin no 
T, ==), Sim PX Cos Oyreos nd 
Te == | COS PKS Time Oyecos oie 
To ==) s€os (Px Cos Oy Sing 
18 = +) Sin PX sim’ Oy Sinene (4 


To save computation time it is necessary to determine 
tren of the eight integrals are always zéro and for whichen 
Ehey are zero. 


Observe that 


cos no = R In? - R,(e)%) 


e 


R, (cos ee, SED ont 


n 
x 
=k, {|—— +j3 — 
ae 
ee xoty 
ae eee 
= (=>) R  (x+jy)’. (SE) 
2 Zz e 
ety 
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Likewise 


sin no = J al I iy)” 
> oe 7 Su). (OE) 


ee yi 
The first term in brackets is always even and may be made 


part of E(x,y). The second term in brackets may be expanded 


by the binomial expansion with the general term 


n-k 


on * Gp” . (75) 


lire tirst integral of (4E) becomes 
I, = [EG y)cos PX cos “Oy R,(xtjy)” dxdy (8E) 


For I, tO Dewmien - Zee , Re (x+jy)” must be even in x and even 
mimes From (7E) this 1s true only of even n; for odd n, Ty 
Will always be zero. 


Examining the other seven integrals yields 


JEG.y) [cos PX cos OY tcos 10 


+ 4 Sin PX Sine) soda iceedy n even 


[E(x.y) [cos Exe sine Oy sone 
| + ‘sin: Px cos7Oy -cocemenlaxrdy J mode 
(9E) 
Spserving that symmetry about the X- and™@y axis 125 necessany 
mor this integral to be non-zero and integrating only when 
ene Integral iS non-zero permits one Co antegrate only over 


tne first quadrant. 
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APPENDIX: 


SUMMATION CHECK 


tine Signal for a poimt target Was sowie in drt meee De 
H(wot) = fN'(r)t(r-e(t)) dér . Cis 
Assuming an infinite image plane, H(w ot) may be written in 
mens O2 the Fourier transform as 


-jamk-p{t) , 
Tak (2F) 


~ 


H(wot) = {N'(k)t*(k)e 


Beaea rectangular detector (ZF) may be written as 


a ey 
fer O (kk sintwkx Sinthky 


Hie) = 7 | 
: > . oe ss 
fou 2 DPM EK Kot Ky Yo) lor) 
et 
x = ¥2 otkx. 
y = v2 onky 


amermoting the integral may be separated inbomtiwoparts 


bY 2) 
hee j— x(pcosw_t - (r_-x_)) 
H(w _t) = = fe eS deci ee teens Q a ne an 
: 420 
San hy 
ae 
2 72 —= £ = - 
1 fy fers y(osinu t - (ry cae 
A y 
(4F) 
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ae _ oe Aa _ a a a | ae ©. ee en Oe >, 


eamce the integration 1s -~ to +~ only even integrands will 











beemon-zero. Rewriting (4F) 
7x 
- > Sin o. a 
ee x O : : 
H(wt) ~ rs e ees oO (2 x[pcosu t - 9)1) dix 
- 5 sin es 
2 ~y ¥20 V2 : _ _ 
-_| e 5 cos ( : y[psinuw t (ry yo) 1) ay. 
=) 
ime general form of the integral is 
i -ué Sin B 
R = R(B) = | e cea ae COs (CWdur Gove) 
O 
imererentiating with respect to B 
dR(B) _ Le wu COS By ecOs  umdu (7F) 
dB tie Y : 
Eicethe solution to (7F) is known [2]; 
_ (c)* | Bec)? 
aR(B) . inraes j | as 
2VT. 
Puce the constant of integration R(0)--— 0, 
» . (re)? ~~ {wte)? 
03) a ae s dv+f e ser G3 2b 
2v¥7 | 0 O 
ne C 
p = Vzc 
and 
Vee 
OF> 7 


ele 





then 














B-C B+C 
2 2 2 2 
R(B) = 5 : [ ede [ 2 Ge 
va J-C/2 Yn J-C/2 
es Bae B+C 
5 fers (a ) a lege ae ( 9 | (10F) 
where 
Z 2 
Eee) moo e* dy. (11F) 
YT O 


Suoestituting the result (10F) into {(5F) 





H(w t) = ee Ere 7 [ocosw pa eet + u W 
O 4 T 0 X O S26 
re ( [ocosegt - x45] - “ 
¥20 
+ Eat v2 [osinw tain oa h 
O O 4 O oe 


| 
¥ 2 bd a 
“Eri, | = 0) OC Seid et eees nee : 2] : COA ys, 
; ( 0 O y “0 ae | 


The mean square voltage of the background may also be 
solved in a like manner. 


From Part V, the mean square voltage was shown to be 


eg = fIsG01° [FGO1? Wyk) a°k : (13F) 


which for a rectangular detector becomes 


Lee ane 
sin-(ak W) 07" *x 
ae x e : 
aa 2a f 2 2s 
B nk oc Oa) ae ; 
x D 


T0 





| pai Pia: os, 
Tice (ae ES el (14F) 
: no ky B°+ (20k,) sb 


Letting 
tes 271K 


the general form of the integral to be solved is 


Lag. Sane 2) oo tt 
Gil) = = a | aa Abe es Ep ( lo) 
O t2/4 0 aie 
Taking the first and second derivatives one has 
00 eae 
ae. _ 2 || Slice) sae Aes alge) 
u ‘1 t ott 
and 
d?R(u) _ 2 [* cos(2yt) .-07t? 
“eee “GE 2 a a (17F) 
du OY - <u 
The solution to the second derivative is known [2]: 
Us ee 
EGE) a ae pe Erfc (ao- 5— 
q 2 20 20 
+ oe! Erfe (ao + y-) | 20) 
where 
Pini ay 
Erfc(z) = &f e% ay (19F) 
Vane 


is the complimentary error function. 
With lengthy and tedious calculations it can be shown 


By USiIn? “tie relation 


oe | 








cr A. 


rp _~ J - CL -— ate u™ © 7 = 


O O 


u 
Ru) =f av L Fwd | (20F) 


tet initial conditions equal to zero and 


2 
E(w) = o Ruy (21F) 
du 


eine t 


Z 
R(v) a i Ente (== - mat 7 (nf2a) 
OL YTU 


le 
ofmmre’ 


S OU u 
+ en c Evie (ao € 76 
Ol 


+e “" Erfc’' (ad - 5) 
= 4 ede (20)] (22F) 


The mean square voltage thus becomes 


et = daBos R(w3a)R(h;8) . et 
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CONTINUE 
SUM = 2.D0:7BF3(1)/X-BF3(2) 
DO 76 T=2eNY¥2.2 


SUM = SUM#2.D04B8F301)3 
CONTINUE 


COMPUTE NORMALIZING CONSTANT 


SNORM = 1.D00/SUM 
DG 7¢@ T=1.NUM 
BF(I)=SNORM?; BF3 (1) 


OCN ORDER BESSEL FUNCTIONS ARE OOD 


WEtKt 2.60.1) BPCIJ=BFCl I=). Dole 
CONTINUE 

BZ = (2.D00#8F3(1)/X-BF3(2))*SNCRM 
RETURN 

END 
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